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A two-parameter extension of the density-scaled double hybrid approach of Sharkas et al. [J. 
Chem. Phys. 134, 064113 (2011)] is presented. It is based on the explicit treatment of a fraction 
of multideterminantal exact exchange. The connection with conventional double hybrids is made 
when neglecting density scaling in the correlation functional as well as second-order corrections to 
the density. In this context, the fraction a c of second-order M0ller-Plesset (MP2) correlation energy 
is not necessarily equal to the square of the fraction a x of Hartree-Fock exchange. More specifically, 
it is shown that a c < a x , a condition that conventional semi-empirical double hybrids actually 
fulfill. In addition, a new procedure for calculating the orbitals, which has a better justification 
than the one routinely used, is proposed. Referred to as Ai variant, the corresponding double 
hybrid approximation has been tested on a small set consisting of H2, N2, Be2, Mg2 and Ar2. 
Three conventional double hybrids (B2-PLYP, B2GP-PLYP and PBEO-DH) have been considered. 
Potential curves obtained with Ai- and regular double hybrids can, in some cases, differ significantly. 
In particular, for the weakly bound dimers, the Ai variants bind systematically more than the regular 
ones, which is an improvement in many but not all cases. Including density scaling in the correlation 
functionals may of course change the results significantly. Moreover, optimized effective potentials 
(OEPs) based on a partially-interacting system could also be used to generate proper orbitals. Work 
is currently in progress in those directions. 



I. INTRODUCTION 

Sharkas et al? recently proposed a rigorous formula- 
tion of one-parameter double hybrid density-functionals, 
which is based on the combination of second-order 
M0ller-Plesset (MP2) perturbation theory with a 
density-scaled correlation functional. The connection 
with regular double hybrids 7 could be made when 
neglecting density scaling. Note that, in their approach, 
the fraction a c of MP2 correlation energy equals the 
square of the fraction a x of Hartree-Fock (HF) exchange. 
Such a condition is actually not fulfilled by conventional 
double hybrids. The a c parameter is either chosen 
irrespective of a x or expressed in terms of a x but not 
as a c — a x . In the latter case, one can for example 
refer to the Perdew-Burke-Ernzerhof-zero double hybrid 
(PBEO-DH) functional of Bremond and Adamo 7 , which 
is characterized by a c 



This work deals with the rigorous formulation of 
density-scaled two-parameter double hybrids. It is or- 
ganized as follows: an exact two-parameter energy ex- 
pression is firstly derived in Sec. |II A| Approximate for- 
mulations are then investigated for defining single hy- 
brid (Sec. II B) and double hybrid (Sec. II C) energy ex- 
pressions. In Sec. II D the connection with conventional 



double hybrids is made. In addition, a new procedure 
for computing orbitals is proposed, defining thus what 
is referred to as Ai variant of the double hybrids. A 
summary of the different approximations that have been 
formulated is then give n in Sec. |II E| Following the com- 
putational details (Sec. Ill), results obtained with regular 



and Ai-double hybrids on a small test set consisting of 
H 2 , N 2 , Be 2 , Mg 2 and Ar 2 are presented and discussed 
in Sec. El 



II. THEORY 
A. Multideterminantal exact exchange 

The approach recently proposed by Sharkas et al. 7 for 
deriving rigorous one-parameter double-hybrid function- 
als is based on the separation of the universal Hohenberg- 

Kohn 7 functional F[n] = F Xl [n] 
partially-interacting contribution 



F Xl [n} = min(*|T + AiWeel*) 



(1) 



with < Ai < 1, and the complement Ai-dependent 
Hartree-exchange-correlation (Hxc) density-functional 



4 1 X cN = (l-Ai)^H X [n]+^ 1 [n], 



E X c 1 [n]=E c {n]-E^H 



(2) 



where -Ehx[«] = [n]|VF ee |<I> KS [n]) is the usual exact 
Hartree-exchange (Hx) term based on the non-interacting 
Kohn-Sham (KS) determinant 7 . The correlation func- 
tionals E c [n] and E^ 1 [n] correspond to fully and partially 



2 



Ai-interacting systems, respectively. The exact ground- 
state energy is then expressed as follows, according to 
the variational principle 7 , 



E = min 



{(* Al NT + X 1 W ee + V nc \^[n}} 



(3) 



where V nc = f dr v ne (r) n(r) is the nuclear potential op- 
erator. Rewritten in terms of a minimization over wave 
functions, Eq. ^ becomes 



E = mm { (*|T + X{W ee + V ae \V) + E^n*} } 



(4) 



= (* A i \f + X X W ee + Ke|* Al ) + £ H xcKm], 

where *Sf Xl fulfills the self-consistent equation 



f + \ 1 W ee + V ne + V Sxc [nyx 1 ] ) |* Al ) = £ Al |* Al ), 



^HxcM 



(5) 



In order to introduce a second scaling factor A2, and thus 
derive two-parameter double hybrids, let us consider the 
following partitioning of the complement Ai-Hxc func- 
tional: 



E Hxc W = ^Hx,mdN + E c!md[n], 



(6) 



which is based on the multideterminantal definition of the 
exact exchange (mdEXX), as introduced by Toulouse, 
Gori-Giorgi and Savin 7 7 7 in the context of range- 
separated density- functional theory (DFT), 

^H X , md W = (1 - X 1 )(^[n]\W ee \^[n]). (7) 

The corresponding complement correlation functional dif- 
fers therefore from the one given in Eq. ([2]): according to 
Eqs. (|6| and Q it can be expressed as 

S c A , 1 md [n]=S c Al [n]-(l-A 1 )A c ^[n], 

A C A >] = ((^[n]\W ee \^[n]) ( 8 ) 
-($ KS [n]|# ee |$ KS [n]) S 



Recombined with Eq. ([3]), Eq. |7| leads to the explicit 
calculation of 100 % of mdEXX. In the spirit of usual hy- 
brid functionals, we want to keep the flexibility of treat- 
ing only a fraction A2 (0 < A2 < 1) of mdEXX explicitly. 
For that purpose, we split the complement mdHx func- 
tional as follows: 



^hW»] = (*a - Xi)(^ Xl [n}\W ee \^[n}) 
+ (1- A 2 )(* A >]|# ee |* Al [n]), 



(9) 



and introduce the complement Ai- and A2-dependent Hxc 
density-functional energy: 



E 



Hxc^W = (1 - \2)(^[n]\W ee \^{n])+E X c : nl M 



= -(X 2 -X 1 ){^[n}\W ee \^[n}) (10) 

+ £ Hx,mdN + E c'md[4 

which, according to Eqs. ^ and ^ , can be rewritten as 



E 



Ai ,A2 



11, A = -(Aa - \i){^[n]\W ee \^[n\) + E%Jn] 

= (l-X 2 )E Hx [n]+E Xl [n] (H) 
+ (X 1 -X 2 )A^[n}. 



Using the following expressions based on a uniform coor- 
dinate scaling of the density 7 7 : 



Et[n} = X 2 1 E c [n 1/Xl ], 
A A 



dE x \ 



ox 



(12) 



A=Ai 



n 1/Al (r) = (l/A 1 ) 3 n(r/A 1 ), 



as well as Eq. (pi), the more explicit density-scaled two- 
parameter (DS2) form 



£hxc A2 M = (1 - A 2 )S Hx [n] + E c [n] - A?£? c [n 1/x J 
+2\ 1 (X 1 - X 2 )E c [n 1/Xl ] 
dE c [n 1/x ] 



(13) 



+Af(Ai-A 2 ) 



9A 



A=Ai 



is obtained and the exact ground-state energy in Eq. ^ 
can then be rewritten as follows: 



E = min 



{(* Al [n]\f+ X 2 W ee + V nc \V Xl [n]> 

;Ai 



-Etc' A} 



(14) 



or, in terms of minimization over local potentials, 
E = mm\ (^ Xl [- ] ^ 1 wJ/ 1 T> l,T,Al 



[v]\T + X 2 W ee + V nc \V M [v]) 

^Ai,A 2 



^Hxc'IVilu]]} > 



(15) 



where x & Al [u] denotes the ground state of the Ai- 
interacting Hamiltonian T + X±W ee + J drv(r) n(r). It is 
important to notice that, for a fixed Ai value, the exact 
minimizing potential v Xl is the one which ensures that 
the density of ^ Al [v Xl ] equals the exact ground-state den- 
sity. Therefore, it does not depend on X 2 . When choosing 
X 2 = Ai, the minimization over potentials in Eq. ( 15 1 can 



be replaced by a minimization over wave functions so that 
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Eq. Q is recovered and, according to Eq. the exact 
potential u Al can be expressed as 



Ai 



oro(r) 



(16) 



As a result, for any Ai and A2 values, the exact ground- 
state energy can be written as 

E=(^\f + X 2 W ee + V nc \^)+E^ 2 [n^ 1 ]. (17) 

If A2 7^ Ai, the energy cannot be obtained straight - 
fowardly from a minimization over wave functions: 



E + min((*|f 



A 2 W e , 

Ai,A 5 



-E 



Hxc 



(18) 



since the minimizing wave function would be eigen- 
function of a A2-interacting system and therefore could 
not be equal to *5> Xl , which is eigenfunction of a Ai- 
interacting system. As discussed in Sec. |IIB| it is in 
principle possible to adapt regular optimized effective 
potential (OEP) methods to this context in order to 
implement Eq. (|15|. 



B. Density-scaled two-parameter single hybrids 

This section deals with the formulation of single hy- 
brid functionals based on the energy expression given in 
Eq. (15). The latter can be considered as a generaliza- 



tion of the standard OEP approach, which is based on 
a non-interacting KS system, since \P Al [v] is the ground 
state of a partially-interacting system. In principle, any 
wave-function-theory-based model can thus be combined 
rigorously with OEPs. We consider in this work the op- 
timization of the potential using the HF approximation 
$ Al [i'] to the exact wave function ^ Al [z;]. Following the 
terminology of Sharkas et al? , we refer to this approach 
as DS2-HF-OEP approximation. In this case, the energy 
becomes 



-&HF-OEP — Him fi H F-OEP |vj 



££E& Aa M = {(* Xi Iv}\t + x 2 w ee + v ne \^[v]) 



(19) 



Such an approach can be interpreted as an exchange- 
only-type OEP calculation based on a Ai-interacting sys- 
tem. As a result, the DS2-HF-OEP energy gradient 
has no simple analytical expression, like in the KS-OEP 
scheme 7 . However, it can be efficiently computed from a 
HF-type linear response vector, using a decomposition of 
the potential in an auxiliary basis of gaussian functions. 



Work is currently in progress in this direction. If we 
simply denote $ and n the converged DS2-HF-OEP de- 
terminant and its density, the energy can be reexpressed, 
according to Eq. (13), as 



E^% X £ X2 = ($|T + V> nc |$) + E B [n] + A 2 £» F [$] 
+(1 - X 2 )E x [n] + E c [n] - Ai(2A 2 - X 1 )E c [n 1/Xl ] (20) 
dE c [n 1/x ] 



+A?(Ai -A 2 ) 



dX 



A=A, 



where E™ [$] is the HF exchange energy. We thus obtain 
a density-scaled two-parameter hybrid (DS2H) exchange- 
correlation energy: 



EC 



= A 2 £« F [$] + (1 - X 2 )EM + E c [n] 

-Ai(2A 2 -A 1 ) J B c [n 1/Al ] (21) 
dE c [n 1/x ] 



+Xf(X 1 -X 2 ) 



dX 



A=Ai 



Note that, in the particular case Ai = A 2 , the exact DS2- 
HF-OEP potential can be expressed in terms of a func- 



tional derivative like in the exact theory (see Sec. II A ) 



Indeed, the DS2-Hxc functional in Eq. ( 13 ) reduces then 



to a DSl-Hxc functional which corresponds to the com- 
plement Ai-Hxc functional: 



£hxc Ai N = (1 - Ai)£Hx[n] + E c [n] - A?£ c [n 1/A J 



(22) 



Moreover, using the one-to-one correspondance between 
the density n$ of any single determinant $ and the local 
potential v Xl [n$] such that the density of the HF deter- 
minant $ Al [v Xl equals n$, the energy in Eq. ( 19 ) can 



be simply obtained by minimization over single determi- 
nants, that is without using OEPs. It then corresponds 
to the DS1H energy of Sharkas et ah' which equals 

£ Al 1H = mm {(*|f + X t W ee + V m \*) + £h*cM} (23) 

The minimizing determinant <E> Al fulfills the self- 
consistent equation 



~ Ai 



r + A 1 tf HF [$ Al ]+Ke + V Hxc [ri^ 1 ] ) |$ Al ) =£*U<^ 1 ), 

(24) 

where C/hf[ ( I >Ai ] is the nonlocal HF potential, so that, in 
this particular case, the minimizing potential in Eq. ( 19 ) 
can be expressed as 



^dsihM = v nc (r) + 



5E 



Hxc 1 



5n(r) 



(25) 
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and the exchange-correlation in Eq. (21) reduces to 



*£Uih = Ai^ F [$] + (1 - A x )£ x [n] + E c [n] 
-X\E c [n 1/Xl ]. 



(26) 



C. Density-scaled two-parameter double hybrids 

This section deals with the formulation of double hy- 
brid fun ctionals based on the energy expression given in 



Eq. (15). In the following, we denote Vq 1 any approxima- 



tion to the exact minimizing potential v Xl . An approxi- 
mate expression for the energy is thus obtained: 

E x^M = + x 2 W ee + V ac \*fr) +^H X c a [n^ 1 ](27) 

where ^q 1 i s the ground state of the Ai-interacting sys- 
tem defined by Vq 



Ai. 



T + XiW e , 



>Ai 



V r 



I* 



= ^o 1 l*o 1 ). 



(28) 



witht^ 1 = J dru^(r)n(r). The energy in Eq. A27l can 
be rewritten as 



T7iAi,A2 c*Ai 
M/n — o n 



(X 2 -X 1 )(^\W ee \^) 
(^iKc-^ 1 !^ 1 ) 



(29) 



which is convenient when applying, as we propose in the 
following, MP perturbation theory to the Ai-interacting 
system described by Eq. ( 28 1 . The zeroth-order wave 



function is chosen to be the determinant $ X which fulfills 
the HF-type equation 



T + X 1 U HF [n 



An 



Ai 



|4> 



An 



= ^|^>. (30) 



we obtain the final energy expression through second or- 
der: 

E* lM = (^\T+X 2 W ee +V ne \$fr) 

+^h 1 X c 2 ["$ a i] + Ai(2A 2 - \i)E^l, 



6E 



(33) 



where Sn^(r) denotes the second-order correction to the 
density associated to ^/ Al . If we simply denote $ and n 
the zeroth-order determinant $ Al and its density, respec- 
tively, the exchange-correlation energy has, according to 
Eq. (13), the form of a density scaled two-parameter dou- 



ble hybrid (DS2DH) functional: 



E 



XC.DS2DH 



X 2 E™[<f>] + (l-X 2 )E x [n] + E [n] 
-Ai(2A 2 -Ai)£ c [ni Al ] 
dE c [m/x\ 



+xj(x 1 -x 2 ) 



+X 1 (2X 2 -X 1 )E- 



dX 



(2) 
MP 



A=Ai 



(34) 



dr 



Une (r)-<(r) + ^^-[n]j^ 2 )(r). 



Let us consider the particular case Ai — X 2 . According to 



Eqs. ( 22 ) and ( 25 ) , when choosing the DS1H determinant 
$ Al and potential v^sm as ^ an d ?V> the second-order- 



density-correction term in the right-hand side of Eq. ( 34 ) 



cancels out and the DS1DH exchange-correlation energy 
of Sharkas et al? is recovered: 



Note that the Brillouin theorem is fulfilled in this con- 
text, which means that the wave function contains only 
double excitations through first order and, as a result, 
the density remains unchanged through first order 7 . The 
perturbation expansion of the fictitious Ai-interacting en- 
ergy through second order equals: 



£ Al = ($*|f + A 1 W ee + V r Al |^ 1 ) + A?£ 



(2) 

where fij^j, is the conventional MP2 energy correction 
calculated with the Ai-interacting orbitals and orbital 
energies obtained from Eq. ( 30 1 . Using the simplified 



2p(2) . 
MP 



(31) 



perturbation expansion for the second term in the right- 
hand side of Eq. p9l: 



(^\W ee \^) = (^\W ee \% 



,Ai\ 



,Ai| 



An 



2Ai-Emp 



(32) 



E M 

-^XC.DSIDH 



A 1 ^ x IF [*] + (l-A 1 )S x [n]+^ c [n] 
-XlE c [n 1/Xl ] + XlE^l 



(35) 



As shown in Sec. II D the DS2DH functional defined in 



Eq. (34) can be connected with conventional double hy- 



brids when neglecting both second-order corrections to 
the density as well as the density scaling in the correla- 
tion functional. 



D. Connection with conventional double hybrids 

In order to connect regular double hybrids with the 
one derived in Sec. |II C[ we neglect both second-order 
corrections to the density as well as the density scaling 
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in the correlation functional: 

<5n( 2 )(r)«0, E c [n 1/X ] « E c [n]. 



(36) 



The DS2H and DS2DH exchange-correlation energies in 
Eqs. (21) and (34 1 reduce then to two-parameter hybrid 



(2H) and two-parameter double hybrid (2DH) exchange- 
correlation energies, respectively: 



and 



£xc x ,2°h = a x E™[<S>] + (1 - a x )E x [n] 
+(1 - a c )E c [n], 

C^h = a x E^[<f] + (1 - a x )E x [n] 
+(1 - a c )E c [n] + a c E^, 



(37) 



(38) 



where we have introduced the two parameters a x and a c 
defined as follows: 



A 2 



Ai(2A a -Ai) 



(39) 



K - Ai) 2 



Let us first notice that, by contrast to KS second-order 
perturbation theory 7 (KS-PT2), single excitations do 
not appear in the 2DH energy expression. Indeed, the lat- 
ter originates from a MP-type calculation, where the Bril- 
louin theorem therefore applies, and which is based on a 
partially Ai-interacting system. In addition, it is readily 
seen from Eq. (39) that, in this context, a c < et x . Inter- 



estingly, the various conventional double hybrids consid- 
ered in Table U fulfill this condition. Note that the one- 
parameter double hybrid (1DH) approximation 7 is recov- 
ered when Ai = A 2 , or, in terms of a x and a c , a c = a x . 
For given values of Ai and A 2 , a unique set of a x and a c 
parameters can be defined according to Eq. ( 39 ) . Recip- 



rocally, for given a x and a c values, we can define the two 
scaling factors Ai and A 2 as follows: 



A 2 



Ai 



(40) 



When a c = 0, the correlation energy is fully described by 
the correlation functional which means that the fictitious 
Ai-interacting system should be the non-interacting (KS) 
one, that is Ai = 0. The standard exchange-only-type 
KS-OEP scheme with a x as fraction of exact exchange is 
thus recovered. We therefore conclude 



Ai = a x - y/aY 



(41) 



Interestingly, the linearly scaled one-parameter double hy- 
brid functional (LS1DH) derived recently by Toulouse et 
al? , which is characterized by a c = a x , is recovered here, 
though density scaling is neglected, when Ai equals 



Af 



a x (l - \/l - a x ). 



Let us now focus on the calculation of the orbitals. As 
pointed out in Sec. \HC\ the DS2DH and thus the 2DH 
exchange-correlation energies are based on HF-type or- 
bitals calculated for a Ai-interacting system. A natural 
choice of orbitals would therefore be the DS2-HF-OEP 
ones. Since density scaling is neglected, we will refer to 
them as Ai-OEP-2H orbitals and the corresponding 2DH 
energy will be referred to as Ai-OEP-2DH. The Ai-OEP- 
2H scheme, which is an OEP-type calculation, can be 
formulated as follows, according to Eqs. (13), (19) and 
(pi and p9l: 



- C/ A 1 -OEP-2H — mill H/ x _OEP-2H [ V \ 1 
V 



(43) 



+E H [n^ 1[v] ] +a x E^[^[v}} 
+(1 - a x )£ , x [n $ A 1[l , ] ] + (1 - a c )i? c [n $ A lM ]|, 
where $ Al [v] fulfills the Ai-interacting HF equation 
[f + \ 1 Ukf[<5> Xi [v]]+v) \^[v])=£^ f [v]\^[v}), 



(44) 

with V = J drv(r)n(r) and Ai defined in Eq. (41). On 
the hand, for conventional double hybrids, the orbitals 
are obtained from a 2H calculation which can be formu- 
lated as 7 



£ 2 a H ,ac = mm{($|f + V> nc |$) + £hM + a*E™m 
+(1 - a x )E x [n$] + (1 - a c )E c [n$] |. 



(45) 



The minimizing determinant $ in Eq. ( 45 ) fulfills the 2H 
equation 



(f + V nc + a x U HF [$] + (1 - ^)Vfex[n 4 ] 



(46) 



+(1 - Oc)Vc[ni] m=S m \9), 



and, therefore, can be interpreted as the HF determinant 
associated to a A 2 -interacting system since, according to 
Eq. (40), A 2 = a x . However, as shown in Sec. II C the 



2DH exchange-correlation functional in Eq. (38) can be 



justified when it is based on HF-type orbitals associated 
to a Ai-interacting Hamiltonian. This condition ensures 
that single excitations do not contribute to the exchange- 
correlation energy. In fact, the conventional calculation 
of the orbitals is simply deduced from Eq. ( 43 ) when re- 



placing the minimization over local potentials by a mini 
mization over single determinants, which is justifie d onl y 
when A 2 = Ai that is equivalent to a c = a x (see Sec 
As shown in Table [TJ conventional double hybrids 



IIBl 



o not 



(42) fulfill the latter condition. In this respect, double hybrids 
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based on the Ai-OEP-2H approximation have a better 
justification than the regular ones. Let us finally con- 
sider a possible alternative to the Ai-OEP-2H scheme, 
which would not require the calculation of an OEP and 
would have a computational cost similar to conventional 
2H calculations. As shown in Sec. |II A[ in the exact the- 
ory and for a fixed Ai value, the OEP should not depend 
on A2. As an approximation, referred to as Ai-2H in 
the following, we assume that this statement still holds 
within the approximate Ai-OEP-2H scheme. As a result, 
an approximate potential can be obtained when choosing, 
in Eq. ( 19 ), A2 = Ai instead of A2 = a x . In this particular 
case, according to Eqs. (1221), plf, and (1361), the Ai-2H 



and 



determinant $' fulfills the modified 2H equation 
(f + K c + a x l> HF [$'] + (1 - ax)t> Hx [n 4 ,] 



+(l-o' c )K[ni,] |*')=3h|3') 



(47) 



where the regular a x and a c parameters have been re- 
placed by Ai and \\ , respectively: 



a' x = Ai = a x — \J a\ — 



a c -> a' c = Af = 2a x fa x - \Ja^- a^j - a c 



(48) 



We thus ensure that the orbitals are computed from a 
Ai-interacting system. This procedure basically consists 
in approximating the exact potential v Xl by 



Wa^hM = w„ c (r) + (1 - 

, , . 5E C 



(49) 



Using this potential and $' as v^ 1 and ^q 1 in Eq. (33 1, we 
obtain from Eqs. (36) and (39) the following expressions 



(50) 



for the Ai-2DH energy: 

E%£S* = (&\? + V nc m + E H [n>] + a x E™[&] 
+(1 - a x )E x [n'] 
+(1 - a c )E c {n'] + a c E^, 



p(2)/ 

where n denotes the density of $ and E\Jp is the regu- 
lar MP2 energy correction calculated with Ai-2H orbitals 
and orbital energies. The Ai-2H and Ai-2DH exchange- 
correlation energies are therefore expressed respectively 



£ X c x ;m%h = a x E™[&] + (1 - a x )E x [n'} 
+{l-a c )E c [h'}, 



(51) 



Ct»H = a x E™[&] + (1 - a x )E x [n'} 
+(l-a c )E c [n']+a c E^. 



(52) 



Let us stress that the Ai-2H and Ai-2DH energies as de- 
fined in Eq. ( 50 ) , and thus the corresponding exchange- 



correlation energies in Eqs. (51) and (52 1, are obtained 



from the same energy expressions as in the regular 2H 
and 2DH schemes, using standard a x and a c values. The 
difference comes from the orbitals which are calculated 
with the modified a' and a' coefficients. 



E. Summary 

A two-parameter extension of the DS1DH scheme pro- 
posed recently by Sharkas et al? has been derived. It is 
based on the explicit treatment of a fraction of mdEXX 
which requires, in the general case, the calculation of an 
OEP for a partially-interacting system. Computing this 
OEP at the HF level of approximation leads to the DS2- 
HF-OEP scheme where the energy has a DS2H form. 
In addition, it was shown that, using any approximate 
potential, DS2DHs can be defined. The connection be- 
tween the latter and regular double hybrids is made when 
neglecting both second-order corrections to the density 
and density scaling. In this case, the DS2-HF-OEP ap- 
proximation reduces to the Ai-OEP-2H one and the cor- 
responding DS2DH scheme (which is based on the Ai- 
OEP-2H potential) is then referred to as Ai-OEP-2DH. 
As an alternative to the Ai-OEP-2H calculation, the Ai- 
2H scheme, where the OEP calculation is replaced by 
a 2H one with modified exchange and correlation coef- 
ficients, has been proposed. A Ai-2DH approximation 
could thus be defined from the Ai-OEP-2DH energy ex- 
pression using Ai-2H orbitals instead of the Ai-OEP-2H 
ones. In this work, only Ai-2H and Ai-2DH results will 
be shown. The OEP-based schemes, which are currently 
under implementation, will be presented in a separate 
paper. 



III. COMPUTATIONAL DETAILS 



Both Ai-2H and Ai-2DH energy expressions in Eq. (50) 



can be simply implemented using a regular DFT code 
that can perform double hybrid calculations. A develop- 
ment version 7 of the DALTON program package 7 has 
been used in this work. In a first step, a Ai-2H calculation 
is done. It consists in performing a 2H calculation with 
the modified exchange and correlation coefficients defined 
in Eq. (48). The corresponding orbitals and orbital ener- 



gies are then used, in a second step, to compute the scaled 
MP2 term. In a third step, the Ai-2H energy is obtained 
from a regular one-iteration 2H calculation, using the 
Ai-2H orbitals as starting orbitals. Three conventional 



7 



double hybrids have been considered: the B2-PLYP 7 
and B2GP-PLYP 7 functionals based on the Becke 88 
(B) exchange functional and the Lee- Yang-Parr (LYP) 
correlation functional, and the PBEO-DH 7 functional 
which is based on the Perdew-Burke-Ernzerhof (PBE) 
exchange functional 7 and the Perdew-Wang (PW) cor- 
relation functional 7 . The corresponding exchange and 
correlation coefficients are given in Table |TJ Calculations 
have been performed on a test set consisting of H2, N2, 
Be2, Mg2 and Ar2. The following basis sets were used: 
cc-pVQZ 7 for H 2 and N 2 , aug-cc-pVQZ 7 for Be 2 and 
Mg2, and aug-cc-pVTZ 7 for Ar2. 



IV. RESULTS AND DISCUSSION 

In this section we compare, in terms of accuracy, the 
regular B2-PLYP, B2GP-PLYP and PBEO-DH double 
hybrids with their Ai variants, as defined in Sec. |IID| 
The potential curves computed for H 2 and N 2 (see Fig.[lJ) 
show a systematic lowering of the total energy when 
using Ai-double hybrids instead of the regular ones. 
This can be analyzed when decomposing the difference 
between the Ai-2DH and regular double hybrid energies 
into a single hybrid term AiSgH = E x *l 2 n — c 
(see Eqs. ((45]) and ( |50| )) and a scaled MP2 term 

a c A£ , MP2 



et al. 



7 study, results are expected to change significantly. 



MP 



E$p) (see Eq. (|50b). As illus- 



trated in Fig. [2] (a), for H 2 , the Ai-2H energy is always 
greater than the regular 2H one (AE'sh > 0) which is 
due to the fact that the Ai-2H orbitals are not optimized 
for the regular 2H energy expression (based on a x and 
a c ) but for a modified one (based on a' x and a' c ). This 
positive difference is then compensated when adding the 
scaled MP2 contribution. According to Eq. (48), the 



fraction a' x of HF exchange used to compute the Ai-2H 
orbitals is lower than the regular fraction a x . This leads 
to a smaller HOMO-LUMO gap (see Fig. [2] (b)) and 
thus to a larger scaled MP2 correction. Note, however, 
that B2GP-PLYP gives the largest et c Ai?MP2 term in 
absolute value while PBEO-DH, for which a x differs the 
most from a x , has the largest reduction (in absolute 
value) of HOMO-LUMO gap. It could be explained 
by the fact that the scaled MP2 correction is sensitive 
not only to the orbital energies but also to the orbitals. 
Results obtained for N 2 (not shown) lead to the same 
conclusion. Returning to the potential curves in Fig. [T] 
we finally note that the Ai variant slightly improves the 
accuracy of the double hybrids in the vicinity of the 
equilibrium distance. 

Concerning the weakly-bound systems Be2, Mg2 and 
Ar2 , the interaction energy curves in Fig. [3] show that 
regular double hybrids and their Ai variants behave 
similarly at large distances but the latter bind more than 
the former. This often leads to slightly more accurate 
results unless the regular double hybrid, like PBEO-DH 
for Be2 and Mg2, already overbinds. It would of course 
be interesting to evaluate the effects of density scaling in 
the correlation functional. Indeed, in the light of Sharkas 



V. CONCLUSIONS 

A rigorous derivation of two-parameter double hybrids 
(2DHs) has been presented. It is based on the combi- 
nation of the DS1DH scheme of Sharkas et al. 7 with 
the explicit treatment of a fraction of multideterminan- 
tal exact exchange. The connection with regular double 
hybrids is made when neglecting both density scaling and 
second-order corrections to the density. It then appears, 
in this context, that the fraction of second-order M0ller- 
Plesset (MP2) energy correlation is smaller or equal to 
the square of the fraction of Hartree Fock (HF) exchange. 
Interestingly, various conventional semi-empirical double 
hybrids fulfill this condition. In the light of those deriva- 
tions, a new procedure for calculating the orbitals, which 
is more justified than the one used routinely, has been 
proposed. It still consists in performing a two-parameter 
hybrid calculation, but with modified exchange and cor- 
relation coefficients. Preliminary results presented in this 
work show that, in such a scheme which is referred to 
as Ai-2DH, the MP2 energy contribution is, in absolute 
value, larger than the regular one. As a result, Ai-2DH 
and regular double hybrid potential curves can, in some 
cases, differ significantly. In particular, for the tested 
weakly bound dimers, the Ai variants bind systemati- 
cally more than the regular ones, which is often but not 
always an improvement. Including density scaling in the 
correlation functionals may of course change the results 
significantly. This still needs to be investigated. More- 
over, optimized effective potentials (OEPs) based on a 
partially-interacting system could also be used to gen- 
erate proper orbitals . Work is currently in progress in 
those directions. 
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FIGURE CAPTIONS 

Figure [TJ Potential curves for H2 (top) and N2 (bot- 
tom). The "exact" curves are taken from Ref. ? 

Figure [2} Single hybrid (A^sh) and scaled MP2 
(a c A£'mp2) energy differences between Ai- and con- 
ventional B2-PLYP, B2GP-PLYP and PBEO-DH 
double hybrids calculated for H2 with respect to the 
bond distance (left); HOMO-LUMO gap for both 
Ai- and conventional B2-PLYP, B2GP-PLYP and 
PBEO-DH double hybrids calculated for H 2 with 
respect to the bond distance (right). 

Figure [3} Interaction energy curves for Be2, Mg2 and 
Ai'2. The accurate 7 , experimental 7 andCCSD(T) 
curves are taken as reference. 
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functional 


d x — A2 


a x 


a c 


Al = fix 


a! c 


B2-PLYP 7 


0.53 


0.28 


0.27 


0.43 


0.19 


B2T-PLYP 7 


0.60 


0.36 


0.31 


0.38 


0.14 


mPW2-PLYP 7 


0.55 


0.30 


0.25 


0.32 


0.10 


mPW2K-PLYP 7 


0.72 


0.52 


0.42 


0.41 


0.17 


B2GP-PLYP 7 


0.65 


0.42 


0.36 


0.40 


0.16 


B2tt-PLYP 7 


0.602 


0.362 


0.273 


0.303 


0.092 


PBEO-DH 7 (o c = o x ) 


0.50 


0.25 


0.125 


0.146 


0.021 



TABLE I: Regular (<z x , a c ) and modified (a x , a' c ) exchange- 
correlation coefficients corresponding to conventional double 
hybrids. 



